On radiative damping in plasma-based accelerators 
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Radiative damping in plasma-based electron accelerators is analyzed. The electron dynamics 
under combined influence of the constant accelerating force and the classical radiation reaction 
force is studied. It is shown that electron acceleration cannot be limited by radiation reaction. 
If initially the accelerating force was stronger than the radiation reaction force then the electron 
acceleration is unlimited. Otherwise the electron is decelerated by radiative damping up to a certain 
instant of time and then accelerated without limits. Regardless of the initial conditions the infinite- 
time asymptotic behavior of an electron is governed by self-similar solution providing unlimited 
acceleration. The relative energy spread induced by the radiative damping decreases with time in 
the infinite-time limit. 
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The plasma-based methods of electron acceleration 
demonstrate an impressive progress in the last ten years. 
The quasimonoenergetic electron bunches are generated 
in laser-plasma acceleration experiments ["J. The elec- 
tron energy in laser wakefield acceleration experiments 
exceeds 1 GeV for cm-scale acceleration length |2| and en- 
ergy doubling of 42 GeV electrons in a meter-scale plasma 
wakefield accelerator is demonstrated Q. Recently the 
physics of linear colliders based on laser-plasma acceler- 
ators have been discussed [4j, |5| . 

The accelerating structure in the plasma-based meth- 
ods is a plasma wave generated behind the driver which 
can be the laser pulse or the electron bunch. There is 
a number of effects which limit the energy gain in the 
plasma-based accelerators [6j. One of the main limita- 
tions comes from the dephasing. The velocity of the rela- 
tivistic electrons becomes slightly higher than the plasma 
wave phase velocity, which is determined by the driver 
velocity. The accelerated electrons slowly outrun the 
plasma wave and leave the accelerating phase. This prob- 
lem can be partially solved by the use of proper longitudi- 
nal gradient of plasma density ~~|, ~§] ■ Another limitation 
is caused by the driver depletion as the driver energy 
converts into the energy of the plasma wave. The driver 
evolution during acceleration (e. g. laser pulse diffraction 
or electron bunch expansion) also imposes certain restric- 
tions on the electron energy gain. In the case of laser- 
plasma accelerators the laser pulse can be guided over 
long distances in the preformed plasma density channel 
["} or with relativistic optical guiding when diffraction is 
compensated by relativistic self- focusing [10]. In general, 
in order to accelerate electrons far beyond the energy 
limited by these effects the multistage schemes can be 
used. 

The electron acceleration in the plasma wave is accom- 
panied with the transverse betatron oscillations caused 
by the action of the focusing force on the electron from 
the plasma wakefield. The accelerating force and the 
focusing force acting on the relativistic electron near the 
driver axis can be approximated as follows F acc — fmcuip 



and Fj_ ~ -nwrii^r, respectively, where r is the trans- 
verse displacement of the electron from the driver axis, / 
and k are the numerical factor and the focusing constant, 
respectively, determined by the parameters of the driver 
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and the plasma, uj p — (4:ire 2 n/niJ is the plasma fre- 
quency, n is the density of the background plasma, m and 
e = — |e| are the electron mass and the electron charge, 
respectively ~§|. For example, if the driver is the lin- 
early polarized Gaussian laser pulse with resonant pulse 
duration then / = 0.35a 2 , ~ 0.7 and n 2 ~ 0.11, where 
a,Q = eEr,/ '(mcwr,) = 2 1 / 2 is chosen, Er is the laser field 
amplitude, wl is the laser frequency ~5|. The period of 
the betatron oscillations is ujp — o; p K7~ 1 / 2 , where 7 is 
the relativistic gamma-factor of the electron. 

The electrons undergoing betatron oscillations emit 



synchrotron radiation 
be estimated as follows P ra d 
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111 Il2j . The radiated power can 
~ 2r e j 2 F 2 /(3mc), where 
cm is the classical elec- 
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tron radius, c is the speed of light. Since the power is 
proportional to the square of the electron energy, the 
radiation losses can stop electron acceleration at some 
threshold value of the electron energy. The threshold 
energy can be estimated by balancing the accelerating 
force and the radiation reaction force, F rr f ~ P ra d/c, so 
that j 2 h ~ f/{en A R 2 p), where Rp — k p r is the normal- 
ized amplitude of betatron oscillations, e = 2r e w p /(3c) 
and k p — uj p /c. The threshold energy is ~ 100 GeV for 
/ = 0.7, n = 10 19 cm- 3 and R p = 1 and k 2 = 0.11. 
Therefore the radiative damping may be a serious limi- 
tation of electron acceleration. 

The electron acceleration in plasma with the radiation 
reaction effect has been studied theoretically [4T.l5l.ll~3", 14 1. 
The radiation reaction has been treated as a perturbation 
[l3j ]. The first-order radiative correction to the energy 
gain of the accelerated electron bunch and the energy 
spread induced by radiation emission have been derived 
for the constant accelerating force. The dependence of 
the electron energy on time has been calculated in the 
plasma channel without the accelerating force and with 
the radiation reaction force [l~3] . Here we study the elec- 
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tron acceleration treating the radiation damping unper- 
turbatively and analyzing the infinite-time limit. 

We start from the relativistic equation for electron mo- 
tion in an electromagnetic field with the radiative reac- 
tion force in Landau-Lifshits form 11511 



du l cr e ik 
7-77- = — F u k 
at e 



2r 2 ■ 

e pi 
rad 

6mc 



(1) 



where F l rad = Ff + F^ + F^, F[ = (e/r e ) (dF ik /dx l ) u k u l , 
pi = _ F ti FklU k ; pi = \ FklU i) (F km u m )u\ F lk is the 
electromagnetic field tensor, u k is the 4-velocity of the 
electron. The first term in Eq. (p} corresponds to the 
Lorentz force and the last term corresponds to the radia- 
tion reaction force. We assume that the ultrarelativistic 
electrons (7 3> 1) are accelerated along a;- axis by the 
force F acc » F±v±/c and undergo betatron oscillations 
driven by the focusing force F± ~ —mK 2 Lu 2 y. Under 
our assumptions, F 3 ^> F 2 and the focusing forces 
make a major contribution to the energy losses through 
radiation. It is convenient to introduce new variables 

P = b,/mc)e 1 /2 / l/2 j Y = T = UptK * / f , 

G = jK 2 f~ 2 . Then Eq. ([I]) can be reduced to the form 



dP 

— = -Y- Y 2 PG, 
dT 



dY 
~dT 



P 



^ = 1-Y 2 G 2 , 
dT 



(2) 
(3) 
(4) 



The obtained equations describe the betatron oscillations 
with the radiative damping. The first term on the right- 
hand side of Eq. ((4]) describes the action of the accelerat- 
ing force, while the second term describes the radiative 
damping. 

When the number of the betatron oscillations is large, 
we can use the averaging method (l6| . To do this let us 
introduce a new variable, S, so that 2S = \U\ = Y 2 + 
P 2 /G = R 2 fe ~ 2(Y 2 ) and U exp (i J GT^dT) = 
Y — iG~ x l 2 P. After averaging over the fast time related 
to the betatron oscillations the averaged equations are 



dS_ 

dT 



_1 S_ 

~2G 
dG _ 

dT ~ 



1 - SG 2 . 



(5) 
(6) 



As G > and S > then dS/dt < and the amplitude 
of the betatron oscillations always decreases with time. 
This means that for arbitrary electron energy the beta- 
tron oscillation amplitude will be small enough at certain 
instance of time to be radiation reaction force less than 
the accelerating force. 

At the absence of the accelerating force (/ = 0), it 
follows from Eqs. ([5]) and (HJ) that SG^ 1 ^ 4 = const and 

/ \-4/5 

7 = 70 f 1 + 5ei?^ j u) p t/16 ) , which is in agreement 




FIG. 1: The phase portrait of the system governed by Eqs. §5§ 
and ©. 



with the solution calculated in Ref. [14j |. where Rp,o — 
Rfj(t = 0). At the absence of the radiation reaction (the 
last terms in RHS of Eqs. ([5]) and © are absent) we get 
G = Gq + T, >/GS = const. The radiation reaction effect 
can be treated as a perturbation. To the first order in the 
radiation reaction force the normalized electron energy 

which is in 



is G = Go 



(2/5) 



1 - (G + T) 5/2 . 

agreement with the result obtained in Ref. 

The system of Eqs. (0 and (HJ) has integral of motion 
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1 - 3SG 2 /2 
S 9 / 4 (SG 2 ) 



3/4 



const. 



(7) 



The electron trajectories in the phase space S — G are the 
integral lines determined by Eq. ([7]). The phase portrait 
of the system governed by Eqs. ([5|) and ([6|) is shown in 
Fig. [1] It is seen from Fig. Q] that if initially the acceler- 
ating force is stronger than the radiation reaction force 
(SG 2 < 1) then the electron energy monotonically in- 
creases with time. Otherwise the electron energy decays 
up to the time instance when F acc — F rr f (that corre- 
sponds to SG 2 — 1) and then it monotonically increases 
with time. It is also seen from Fig.[T]that all electron tra- 
jectories merge in the the limit t — ¥ 00 so that G — > 00 
and S -> 0. It follows from Eq. © that S = 2G-' 2 /3 in 
this limit. We will call the electron acceleration in this 
limit as an asymptotic acceleration regime (AAR). 

We verify our analytical results by numerical simula- 
tions. The exact equation ([T]) and the averaged equations 
of motions ([5]) and ([B]) are integrated numerically for test 
electrons for / = 0.1 and n = 10 15 cm -3 . For simplicity, 
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FIG. 2: The dependence of a) 7 and b) I„ on u p t calculated by 
solving of the exact Eq. (fTJ) (black solid lines) and by solving 
of the approximate Eqs. J!])-© (red dashed lines) for / = 0.1, 
k 2 = 0.5, n — 10 15 cm -3 and for initial conditions 70 = 2000, 
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0.8, p y ,o = 0. 



we consider the structure of the transverse electromag- 
netic field similar to the bubble regime: k 2 = 0.5 and 
E±_ ps H± . The dependence of the normalized integral of 
motion /„ = 7 _1 (eK/ 2 ) -3 , and 7 on uj p t for initial con- 
dition 70 = 2000 and Rp^ = 0.8, p y $ = is shown on 
Fig. O It is seen from Fig. [5] that the solution of the ex- 
act equations and that of the approximate averaged are 
in a good agreement. Moreover, the integral I is almost 
constant for the exact equations (pj (see Fig. [2}:). 



We can introduce new variables g — G/Gt r , t = T/T tr 
and s = (S/Str) (G/G ir )~ 1/4 , where G tr = T tr = S^ 2 = 
I 2 / 9 . Then Eqs. ([3]), (JB]) and © are reduced to the form 
which does not depend on any parameters. Therefore 
the characteristic time of transition to AAR is ~ T tr - 
The solution of the equations can be written in term 
of hypergeometric function, 2-P1 (fli b; c; z), [18] as fol- 
lows <p(a)-<p (s ) =r,ip (s) = 2 4/9 (3 + 2s 2 ) 5/9 s" 4 / 9 - 
2 i3/9 3 5/9 s i4/9 aJ?1 (7/9^/9; 16/9; -2s 2 /3), where s = 
s (r = 0). The asymptotic expansions of function <p(s) 
are ip(s) pa 3(3s/2)~ 4/9 for s < 1, cp(s) ps 8 + s" 4 / 3 
for s 3> 1, where 8 ps 1.85. Thus in the limit t ^> 1 
s ~ r _9/4 < 1 and g ~ r > 1. 

To derive the asymptotic solution the initial condi- 



tion should be applied. We assume that SqGq <C 1 (so 
that s <C 1 and I ~ S~ 3 G~ 3 ^ 2 ) which is typical for 
the initial parameters of the electron beam. For exam- 
ple, this condition is fulfilled for the initial parameters 
- 2 < 0.1 TeV, n < 10 18 cm" 3 , Rp >0 = 1, / = 0.7, 
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K = 0.11. Making of use the asymptotic expansion for 
s < 1 and s > 1 we have (9/4)s _4/9 ps t + 8. There- 
fore the normalized electron energy and the square of the 
normalized betatron amplitude are in the limit T ^> T tr 



G 



-Gtr 



-T, S 



-G~ 



(8) 



We can conclude that in AAR F rr f = 2F acc /3 so that 
the electron energy increases linearly with time while the 
betatron amplitude is reversely proportional to the time. 

The averaged equations of motions J5|) and (JB]) are in- 
tegrated numerically for the test electrons with the same 
parameters as for Fig. [5] for three values of the initial 
betatron amplitude Rp,o = 0.8, 0.2, 0.1. It is seen from 
Fig.[3]that the asymptotic solution ([§]) is in a good agree- 
ment with the result of numerical integration. 

The radiation damping rate varies for the electrons 
with different betatron oscillation amplitudes. This 
causes the energy spread in the electron bunch acceler- 
ated in the plasma wave. We assume that the ampli- 
tude of the betatron oscillations of the electrons in the 
accelerated bunch is uniformly distributed in the range 
R min < Rp fi < R m ax and R max > Rmin- We also again 
assume that So G 2 -C 1. Then the normalized mean en- 
ergy and the normalized square of the relative energy 
spread are in AAR 



(G> 



R 2 

'-in 



f T 

/ GRpfidRpfi — G max 8 + — , (9) 



R m ir, 



a 2 G = (G 2 ) - (G) 2 ~ G 2 max — y 



Rn 
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(10) 



where G max = G tr (Ri3fi = Rmax)- It follows from 
Eqs. © and (fTU)) that the relative energy spread, 
ac I (G) , decreases with time in AAR. 

Eqs. ©-((ll) are derived under conditions that F± gives 
the main contribution to the radiative damping and 
F3 Fi, F2. However F± goes to zero in the limit 
t 00. Therefore we should check: should the acceler- 
ating force and terms F± , F% be taken into account in the 
radiation reaction force in this limit? First it is signif- 
icant that the radiation reaction force remains constant 
in AAR because F± ~ Rp ~ > and 7 — > 00 for t — » 00 
in such way that -R«7 2 = const. Making of use Eq. ([8]) 
and relation v y ~ uipy we get F2/F3 ~ ft <C 1 and 
Fi/F 3 ~ (3/4)K 2 / 7 - 1/2 e 1/2 < 1, where we assume that 
K ~ / ~ 1. The contribution from the accelerating force 
(or from E x ) to F3 is of the order F2/F3 <C 1. Therefore 
our model defined by Eqs. (O and ([BJ is valid in AAR. For 
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FIG. 3: The dependence of 7 on u p t in AAR: analytic solution 
(red dashed lines) and numerical solution (black solid line) for 
Rp,o = 0.8 (lines 1), R P:0 = 0.2 (lines 2) and R 0:O = 0.1 (lines 
3). The other parameters are the same as in Fig. [2] 



high energy electrons quantum electrodynamics (QED) 
effects can be important. The energy of the photon emit- 
ted by the accelerated electron can be so high that the 
quantum recoil becomes strong. The photon emission 
can be treated in classical approach if QED parame- 



1/2 



ter x = (mc7E + p x H)" - (p • E)" / (mcE cr ) ~ 
r yF±/(eE cr ) is much less than unity, where E cr — 
m 2 c 3 /(eh) « 1.32 x 10 16 V/cm is the QED critical 
field [l7|. x can be estimated in AAR as follows x ~ 

[(2f/a) (hcu p /mc 2 )] 1/2 < 1, where a = e 2 /hc pa 1/137 
is the fine structure constant. Therefore the classical ap- 
proach for the radiation reaction force is valid in the limit 
t — > 00 because, like for the corrections to the radiation 
reaction force, the growth of 7 in x is compensated by 
decreasing of F±. 

The distance passed by the electron before reaching 

AAR is k p l tr ~ (f/K 2 )T tr ~ 1.6 (e 2 lQ Rl Q fK 8 y . 

For the initial parameters n = 10 18 cm" 3 , Rp = 1, 
70 = 2 • 10 3 , / = 0.7, k 2 = 0.11 the electron comes into 
AAR after passing 7800 laser-driven acceleration stages 
with total distance lt r — 73 m, achieving the energy 
7771c 2 ~ 5 TeV and R@ ~ 0.008, where the stage dis- 
tance is chosen to be equal to the half dephasing length 
H and the distance between the acceleration stages is 
neglected. For the rarefied plasma 77 = 10 15 cm -3 , AAR 
is achieved in 78 stages with l tr ~ 23 km, 7777c 2 ~ 48 TeV 
and Rp ~ 0.005. AAR may be achieved within one accel- 
eration stage in the proton-driven acceleration schemes 
because of very large dephasing length 

In conclusions, we have shown that the electron accel- 
eration is not limited by the radiative damping in plasma- 
based accelerators. Even if the radiation reaction force 
is stronger than the accelerating force at the beginning, 
then acceleration eventually succeeds deceleration with 
time. The damping of the betatron oscillations leads to 



the transition to the self-similar asymptotic acceleration 
regime in the infinite-time limit when the radiation re- 
action force becomes equal to two thirds of the acceler- 
ating force. The relative energy spread induced by the 
radiative damping in the accelerated electron bunch de- 
creases with time in this regime. This opens possibility 
to use high density plasma at the late stages of multi- 
stage plasma-based accelerators despite the fact that the 
radiative damping is enhanced as density increases. The 
high density plasma can be favorable because it provides 
high accelerating gradient and, thus, reduces the length 
of the acceleration stages. The obtained results can be 
also applied to any other accelerating systems with the 
linear focusing forces. 
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